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SUMMARY 


Instead  of  solving  the  nonlinear  differential  equation  which  governs 
the  compressible  flow,  an  approximate  method  of  solution  by  means  of  the 
variational  method  is  used  in  this  report.  The  general  problem  of  steady 
irrotational  flow  past  an  arbitrary  body  is  formulated.  Two  examples 
were  carried  out,  namely,  the  flow  past  a  circular  cylinder  and  the  flow 
past  a  thin  curved  surface.  The  variational  method  yields  results  of 
velocity  and  pressure  distributions  which  compare  excellently  with  those 
found  by  existing  methods.  These  results  indicate  that  the  variational 
method  will  yield  good  approximate  solution  for  flow  past  both  thick  and 
thin  bodies,  at  both  high  and  low  Mach  numbers. 


INTRODUCTION 


The  advance  in  aerodynamics  has  been  greatly  accelerated  by  the 
interpretation  of  good  experimental  results  in  the  light  of  theoretical 
predictions.  When  the  flow  velocity  is  small,  the  change  in  density  of 
the  fluid  in  the  flow  is  so  small  that  the  fluid  may  be  considered  to  be 
incompressible.  Many  useful  theoretical  results  have  been  obtained  from 
a  consideration  of  the  two-dimensional  irrotational  flow  of  an  incom¬ 
pressible  perfect  fluid.  When  the  flow  velocity  becomes  large,  the 
change  in  density  of  the  fluid  in  the  flow  can  no  longer  be  considered 
small  and  the  flow  must  be  considered  to  be  compressible.  While  the  two- 
dimensional  incompressible  flow  can  be  described  by  the  well-known 
Laplace  equation,  the  corresponding  compressible  flow  can  only  be 
described  by  a  complicated  nonlinear  differential  equation.  Except  for 
a  few  simple  cases,  the  exact  solution  of  this  differential  equation  is 
still  unknown. 

In  the  case  of  two-dimensional  irrotational  subsonic  flow,  there 
are  at  least  three  different  methods  of  small  perturbation  or  lineariza¬ 
tion.  The  small -perturbation  methods  have  been  carried  out  by  developing 
the  velocity  potential  in  terms  of  the  Mach  number  or  of  a  thickness 
parameter.  Obviously,  the  development  in  terms  of  the  Mach  number  limits 
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the  applicability  of  the  method  to  flows  at  small  velocities  and  the 
development  in  terms  of  the  thickness  parameter  limits  the  applicability 
of  the  method  only  to  flows  about  very  thin  bodies.  Methods  for  calcu¬ 
lating  the  second-  and  the  third -order  approximations  in  these  methods 
have  been  carried  out  in  some  cases,  -but  the  increased  mathematical 
difficulty  with  each  approximation  prevents  their  application  to  the 
higher  orders  of  approximation.  Another  method  of  linearization  is  the 
so-called  hodograph  method  wherein  the  nonlinear  differential  equation 
is  transformed  into  an  exact  linear  differential  equation,  thus  removing 
the  difficulty  of  the  nonlinearity  of  the  governing  equation.  However, 
in  doing  so,  the  fulfillment  of  the  boundary  conditions  is  made  much 
more  difficult.  The  application  of  this  method  to  actual  problems 
requires  some  simplifying  assumptions  which  can  be  shown  as  being  only 
approximately  correct. 

The  variational  method,  as  carried  out  in  this  report,  on  the  other 
hand,  can  be  used  to  study  the  compressible  flow  without  linearization 
and  the  computation  can  be  made  to  obtain  a  solution  of  high  accuracy 
without  any  mathematical  difficulty.  The  variational  method  has  also 
the  advantage  of  being  a  direct  method  by  means  of  which  one  may  be  able 
to  calculate  the  flow  passing  a  given  airfoil  at  different  Mach  numbers. 
In  the  variational  method,  instead  of  solving  the  nonlinear  differential 
equations  of  motion,  a  variational  integral  is  first  formulated  and 
approximate  solutions  are  then  calculated  by  the  Rayleigh-Ritz  procedure. 

Two  numerical  examples  are  carried  out  in  this  report  to  ascertain 
the  applicability  of  this  method.  They  are  for  the  two-dimensional 
irrotational  flows  past  a  circular  cylinder  and  a  thin  bump  as  proposed 
by  Kaplan.  (See  reference  1.)  These  examples  are  chosen  because  in  the 
case  of  circular  cylinders  similar  results  by  the  Rayleigh-Janzen  method 
are  available,  and  in  the  case  of  the  Kaplan  bump  results  by  means  of 
the  perturbation  method  in  terms  of  thickness  parameter  are  available. 

A  comparison  of  the  results  by  the  present  method  with  those  computed 
by  the  other  methods  shows  excellent  agreement.  Since  the  Rayleigh- 
Janzen  method  is  known  to  give  good  results  for  thick  bodies  at  a  rela¬ 
tively  low  Mach  number  and  the  perturbation  method  in  terms  of  thickness 
parameter  yields  accurate  results  for  thin  bodies  up  to  high  Mach  numbers 
the  results  of  this  report  indicate  that  the  variational  method  gives 
good  approximate  solutions  for  both  thick  and  thin  bodies  and  at  both 
high  and  low  Mach  numbers. 

This  investigation  was  carried  out  at  the  Daniel  Guggenheim  School 
of  Aeronautics,  College  of  Engineering,  New  York  University,  under  the 
sponsorship  and  with  the  financial  assistance  of  the  National  Advisory 
Committee  for  Aeronautics.  The  author  wishes  to  acknowledge  his  indebt¬ 
edness  to  Messrs.  J.  R.  Knudsen  and  H.  H.  Hilton  who  carried  out  most  of 
the  numerical  calculation.  Grateful  acknowledgment  is  hereby  also  made 
to  Professor  K.  0.  Friedrichs  for  his  kind  encouragement  and  inspiring 
discussions  and  to  Professor  F.  K.  Teichmann  for  his  understanding 
cooperation. 
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SYMBOLS 

a  velocity  of  sound 

Cp  pressure  coefficient 

C_  4  pressure  coefficient  for  incompressible  flow 

Pj-1- 

I  variational  integral 

M  Mach  number  (q/a) 

n  _  inward  normal  direction 

p  pressure 

q  velocity  vector 

q  magnitude  of  velocity  vector 

qmax  maximum  possible  velocity  of  flow 

r  position  radius  in  flow  field 

u  velocity  component  in  x-direction 

U  velocity  of  undisturbed  stream 

v  velocity  component  in  y-direction 

x,y  rectangular  coordinates  of  a  point  in  the  fluid 

X,Y  coordinates  of  curved  surface  divided  by  semichord  of  shape 

z  =  x  +  iy  =  rei@ 

t  thickness  parameter 

6  constant  defining  thickness  of  Kaplan*  s  bump 

7  ratio  of  specific  heats  at  constant  pressure  and  constant 

volume 

6  angular  position 

!  velocity  potential 
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p  mass  density  of  fluid 

£  ,t]  rectangular  coordinates  in  plane  of  curved  surface 

C  =  £  +  Itj 

T  strength  of  circulation 

Subscripts: 

x,y  differentiation  in  corresponding  direction 

o  conditions  of  undisturbed  stream 


THEORY 


For  steady  two-dimensional  irrotational  flows,  there  are  four  hydro 
dynamical  variables:  The  pressure  p,  the  density  p,  and  the  two 
components  of  velocity  u  and  v.  All  these  are  to  be  determined  as 
functions  of  the  coordinates  x  and  y  by  the  following  four 
equations : 


(1)  The  equation  of  continuity  (one  equation) 

div  ( pTj)  -  0 


or 


|p  +  =  ° 

where  q  is  the  velocity  vector. 

(2)  The  equations  of  motion  (two  equations) 


Su  ,  Su  -1  S' 

use  +  v^=-p 


Sv  3v  -1 
u  +  v  ^  =  -p 


(3)  The  equation  of  state  (one  equation) 

p  =  A  +  Bpy 


(1) 


(2) 


(3) 


where  A,  B,  and  y  are  constants  and  the  significance  of  A  and  B 
will  be  discussed  later. 
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Most  aerodynamic  problems  are  concerned  with  either  flow  from  rest 
or  flow  which  is  parallel  and  uniform  at  infinity.  Both  of  these  types 
are  irrotational  to  begin  with,  and  will  remain  so  provided  that  fric¬ 
tion  and  shock  waves  are  neglected.  Thus,  irrotational  flow  is  of  great 
interest  in  aeronautical  applications. 

The  condition  of  irrotationality  is 

rot  q  =  0 
or 

<3u  bv  _  n 

5^  ~  S  ” 

This  condition  can  be  satisfied  by  introducing  a  velocity  potential  0 
such  that 


w 


that  is,  q  is  a  gradient  of  0. 

For  the  velocity  field  given  by  equations  (4)  to  represent  a  physi¬ 
cally  possible  field  of  flow^  it  is  necessary  that  the  continuity  equa¬ 
tion  (1)  be  also  satisfied.  By  means  of  equations  (l),  (2),  and  (3),  p 
and  p  are  eliminated,  and  one  obtains  the  following  fundamental  dif¬ 
ferential  equation: 

(a  -  0X  0y  J  0yy  *  20x$y0xy  =  0  (  5) 

in  which 

„2_^2,7-l/  2  ^  2  *2\ 

a  =  ao  +  -  h  ) 

where  a  is  the  velocity  of  sound  in  the  fluid  and  the  subscript  o 
refers  to  quantities  at  infinity.  Equation  (5)  is  a  nonlinear  differ¬ 
ential  equation,  and  its  exact  solution  is  found  to  be  extremely  diffi¬ 
cult.  Instead  of  solving  the  differential  equation  exactly,  an  approxi¬ 
mate  solution  may  be  obtained  by  means  of  the  variational  method.  In 
this  method  the  first  step  is  to  formulate  a  variational  principle. 
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A  variational  principle  for  the  compressible  fluid  was  first 
formulated  by  Hargreaves.  (See  reference  2.)  A  somewhat  more  detailed 
discussion  of  this  variational  principle  was  made  later  by  Bateman. 

(See  reference  3.)  A  simplified  proof  of  the  principle  was  given  by 
Wang.  (See  reference  4.)  It  seems,  however,  that  the  following  formu¬ 
lation  is  more  complete  from  the  mathematical  point  of  view  and  is 
therefore  included  here. 


Consider  the  following  variational  integral 


1  =  /  <f(p) 

°  D 


U  ^  +  V 

ox  <3y 


|(u2  +  V2) 


dx  dy 


where  the  integration  is  taken  over  the  given  domain  D.  Varying  p, 
u }  v,  and  0  independently  and  setting  the  first  variation  equal  to 
zero }  one  obtains 

&1  =  J  *1  f •  (p)  -  (u  H  + v  +  |(u2  +  -  p(H  -  u)5u  - 

p(|^  ’  V)5V  +  (lx*’  +  ly^)5^l  dy  +  J  p8^n  dS 


where  S  is  the  boundary  of  the  domain  D  and  qn  is  the  velocity 
component  which  is  normal  to  the  boundary  S. 

The  Euler's  equations  of  the  integral  I,  according  to  the  calculus 
of  variations,  are  then 


f(p)  4i*vi§ 


+  i(u2  +  v2)  =  0 


-  u  =  0 


v  =  0 


Spu  dpv  _  o 

ST  ST 


(10) 
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and  the  natural  boundary  condition  must  be  such  that 


p&0qn  dS  =  0 


(ID 


Equation  (10)  is  the  continuity  equation  (l)  and  equations  (8)  and  (9) 
are  the  same  as  equations  (4),  which  implies  that  the  condition  of  irro- 
tationality  is  satisfied.  If  the  pressure  p  is  given  by  the  equation 

p  =  -pf'(p)  +  f(p)  (12) 

equation  (7),  when  combined  with  equations  (8),  (9),  and  (10),  indicates 
that  equations  (2)  are  satisfied. 

Thus  the  variational  principle  (6)  leads  to  the  fundamental  equa¬ 
tions  of  a  two-dimensional  compressible  fluid  if  p  is  of  the  form 
given  by  equation  (12)  and  the  condition  (11)  is  satisfied.  Equa¬ 
tion  (12)  is  a  differential  equation  of  the  Clairaut  type  and  its  solu¬ 
tion,  when  p  =  A  +  Bp7,  is 

f(p)  =  A  +  <W2  p/2  "  Bp7A?  ■  (13) 


Where  Vx2  =  7“TT  aQ2  +  %  &2  =  =  ^P7’1-  The  condi¬ 

tion  (ll)  is  satisfied  if  the  domain  is  finite  because  the  boundary  con¬ 
ditions  in  fluid  dynamics  are  either  that  the  boundary  is  a  stream  sur¬ 
face,  thus  qn  =  0,  or  that  the  velocity  is  prescribed;  that  is,  80  =  0. 
For  the  domain  which  extends  to  infinity,  since  a  part  of  the  boundary  S 
is  infinite,  the  integral  may  result  in  a  finite  value,  although  the 
integrand  may  approach  zero  at  infinity.  This  case  has  already  been  dis¬ 
cussed  in  detail  by  Wang  (reference  4)  and  in  such  cases  this  finite 
value  of  integral  (ll)  should  be  subtracted  from  integral  (6). 

In  applying  the  Rayleigh -Ritz  method  to  the  problem,  it  is  more 
convenient  to  use  only  one  variable  in  the  variational  integral,  say  the 
velocity  potential  0,  and  to  express  the  other  variables  u,  v, 
and  p  in  terms  of  0  by  means  of  equations  (8),  (9),  and  (10).  In  so 
doing,  equations  (8),  (9),  and  (10)  are  automatically  satisfied  and  the 
first  variation  of  equation  (6)  then  leads  to  the  fundamental  differential 
equation  (5).  With  equations  (8),  (9),  (10),  and  (12)  satisfied,  it  can 
be  easily  verified  that  the  integrand  of  I  is  just  the  pressure  p. 

The  Rayleigh-Ritz  method  consists  essentially  of  choosing  a  series 
of  functions  which  satisfies  the  boundary  conditions  but  has  undetermined 
parameters,  such  as  A^,  i,j  =  1,  2,  .  .  .,  in  the  functions.  If  the 
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variational  integral  is  denoted  by  I,  then  by  substituting  these 
functions  into  .  I  and  setting  the  first  variation  of  I  equal  to  zero, 
the  following  set  of  simultaneous  algebraic  equations  is  obtained: 

3a  =  0  =  1,  2,  .  •  . 

dAij 

The  parameters  A.^  can  thus  be  determined  by  solving  these  simultaneous 

equations.  The  Rayleigh -Ritz  procedure,  therefore,  gives  a  solution 
which  satisfies  the  boundary  conditions  exactly  but  satisfies  the  dif¬ 
ferential  equations  only  approximately. 

Consider  the  case  of  a  compressible -fluid  flow  past  a  circular 
cylinder  with  unit  radius.  The  velocity  potential  0  in  this  case  may 
be  assumed  to  be 

0  =  0-l  +  02  (14) 

with 

01  =  u(r  +  i)cos  e  +  £  0  (15) 

and 

02  '  £  tl  '  <-  *  2) r**)  CO°  ”9  +  ^  ”*)  <16) 

where  r  is  the  strength  of  circulation,  and  are  undetermined 

parameters,  and  r  and  6  are  the  polar  coordinates.  It  is  evident 
that  for  the  assumed  value  of  0  the  boundary  conditions  are  satisfied; 
that  is,  at  r  =.00  and  q  =  U,  U  being  the  velocity  of  the  undisturbed 
flow  at  infinity  and  on  the  surface  of  the  body,  S0/5 n  =  0  or  the 
normal  component  of  the  velocity  is  zero,  n  being  the  inward  normal 
direction. 

The  circular  cylinder  in  the  complex  z-plane  may  be  mapped  into 
cylindrical  bodies  of  arbitrary  shape  in  the  complex  £ -plane  by  the 
mapping  function  £  =  f(z).  When  the  mapping  is  conformal,  it  was  proved 
in  reference  4  that  0(£)  still  satisfies  the  boundary  conditions. 

With  0  as  assumed  in  equations  (15)  and  (l6),  the  final  variational 
principle  (reference  4)  may  be  written  as  follows: 
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7]7/(7-l) 

I  =  IT  q 1  2 * * * * 7  -  SilL  ^  r  dr  d0  + 

Jj  max  d£  2  dz 


z -plane 


1  ISm 


2  „2 


1/ (7-1) 


where  the  last  term  Is  due  to  the  fact  that  the  integral  (ll)  is  not 
zero  and 


{  ,2  (b$X 

q(z)  =  {&) 


The  final  velocity  in  the  plane  of  this  arbitrary  body  is 

q(02  =  q(z)2  ||  (19 

The  integral  equation  (17)  can  render  itself  into  analytical  inte¬ 
gration  only  when  7/(7  -  1)  is  an  integer.  For  different  values  of 
7/(7  -  1),  the  corresponding  values' of  7  are  as  follows: 

7/(7  -  1)  7 


1  1.333 

The  isentropic  value  of  7  for  air  is  1.405,  which  is  between  the 
values  of  1.5  and  1.333-  Among  these  values  of  7,  7  =  1.333  gives 

the  closest  approximation  to  the  isentropic  value.  However,  in  this 

case  the  integrand  contains  an  expression  which  is  raised  to  the  fourth 
power  and  the  labor  in  carrying  out  the  integration  becomes  excessive. 

In  order  to  simplify  the  numerical  work,  7=2  is  used  in  the  subse¬ 

quent  calculation.  Chaplygin  and  Von  Karma n  and  Tsien  used  7  =  -1  in 
their  solution  by  the  hodograph  method.  The  justification  of  taking 

7  =  -1  has  also  been  discussed  by  the  latter  two  authors.  The  idea  is 
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that  the  velocity  in  the  flow  field  does  not  differ  too  much  from  the 
undisturbed  velocity. and  therefore  only  a  small  portion  of  the  p  p 
curve  will  he  used.  If  one  makes  the  curve  with  7  =  -1  tangent  to  the 
curve  with  7  =  1.405  at  the  point  that  represents  the  conditions  at 
infinity  (pQ,pQ)  and  uses  the  curve  with  7  =  -1,  the  error  will  not  he 
too  large.  The  use  of  7  =  2  instead  of  7  =  -1  gives  even  better 
approximation  as  can  he  readily  seen  from  figure  1.  Although  the  curves 
with  7  =  -1  and  7  =  1.405  may  he  made  tangent  to  each  other  at 
(p  ,pQ),  their  second  derivatives  have  opposite  signs.  When  the  curves 
wi§h  7=2  and  7  =  l.k05  are  made  tangent  to  each  other  at  (p0jPo)> 
however,  a  small  portion  of  the  curves  is  very  close  together,  and  thus 
a  much  better  approximation  is  provided. 

The  isentropic  pres sure -density  relationship  is  as  follows: 

P'  =  po'(p,)1-405/(po,)1*1+05  (20) 

According  to  the  previous  discussion,  the  constants  A  and  B  in 
equation  (3)  can  he  determined  from  the  condition  that  the  curve  given 
by  equation  (20)  and  that  by  equation  (3)  are  tangent  to  each  other  at 
(p0,p0).  This  yields  two  conditions)  namely, 
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For  flow  past  the  circular  cylinder,  the  range  of  variation  of  p/p 
is  approximately  from  0.2  to  1.2.  In  this  region,  it  can  be  easily* 
seen  from  figure  1  that  equation  (23)  does  indeed  give  a  very  good 
approximation  of  equation  (20). 


FLOW  PAST  A  CIRCULAR  CYLINDER 


The  problem  of  determining  the  flow  past  a  circular  cylinder  with¬ 
out  circulation  has  been  carried  out  by  Wang  and  reported  in  reference  4 
A  detailed  discussion  has  been  made  there  for  the  flow  at  mq  =  0.4 
Because  of  the  symmetry  of  the  problem,  0  may  be  assumed  as  follows: 


0  U\v  +  rjcos  9  +  A±1(l  -  ^j-jcos  6  +  Al3^p  -  ^jjcos  36  + 


Mp  ■  i?. rs  6  *  ■  s)cob  38  + 


M?  •  i?)cos  58  *  *51(^5  - 

In  expression  (25),  six  unknown  parameters  are  taken.  The  variational 
integral  is 


pr=oo  p  @= 

r=l  J  6=0 


2  -  ®2  -  (?  S)2 


r  dr  &Q  +  k( q _ 2  -  U2Wa. 


where  7  is  taken  as  2. 

Substituting  0  as  given  by  equation  (25)  into  equation  (26)  and 
differentiating  the  resulting  integral  with  respect  to  the  parameters, 
one  obtains 


3i/6a11  =  o 

blJbA^  =  0 

Al/M31 

A. 

Sly/cJA^,  =  0 

bl^A^  =  0 

aiKi 
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Integrating  these  equations,  there  results 

(2.07^07qmax2  -  7.31852U2)ai;l  +  7.o6667U2a13  +  (0.325926^^  - 
i.0370to2)A31  +  o.915555U2a33  +  (o.ll4286qmax2  -  0.31+0136u2)a51  =  K± 

7.o6667U2a11  +  (13 • 3333qmax2  -  44.2667U2)a13  +  i.20000U2a31  + 

^1 . 60000 qjjjjQj2  -  5.9733^u2)a33  +  32. WrtJ2A15  +  o.44o8i6u2a51  =  k2 

(0.32992 6^^  -  i.03704u2)ai;l  +  i.20000U2a13  +  (0.09007^1^/  - 

0.207407U2 )a31  +  0.19225to2A33  +  (o.0^1l428qmax2  -  0.0810885U2)a51  =  Kj 

■ 

0.915555U2A11  +  (l. 60000 qj^2  -  5-9T33to2)A13  +  0.192254u2a31  + 
(o.277333qmax2  -  1 . 01181U2 )a33  +  4.924WrtJ2A15  +  0.08l3062U2A51  ==  K4 

32.1(.4^U2A13  +  4.92UAu2A33  +  (35.8519qmax2  -  120.9^802)a15  =  Kj 

(0. nit-286^ y2  -  0 . 340136u2)ai;l  +  O.A4o8l6u2A13  +  (o.oAll428qTnay2  - 
0.0810885U2)a31  +  0.0813062U2A33  +  (o.02l8309qmax2  - 
0.0360739U2)a51  =  Kg 


where  Kp,  K2,  K3,  K^,  K^,  and  Kg  are  sums  of  nonlinear  terms  in 

AU,  A13,  A31,  A33,  A15,  and  A^  and  their  complete  expressions 

are  given  in  reference  k,  except  that  the  first  term  in  Kp  should 
be.  4.22222U3  and  the  first  term  in'  K3  should  be  -6.66666u3;  the 

corresponding  values  given  in  reference  k  are  due  to  misprints.  These 
equations  .can  be  solved  as  follows:  First,  assume  the  values  of  the  A's  a 

to  be  zero  in  the  K's.  Equation  (27)  then  becomes  a  system  of  linear 
simultaneous  equations  which  can  be  solved  by  the  method  of  solving 
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linear  simultaneous  equations  suggested  by  Crout  (reference  5) .  The 
values  of  A -y_,  A-^,  A^,  A33,  A51  so  computed  may  be 

taken  as  the  first  approximation.  Use  the  values  of  the  ATs  so  deter¬ 
mined  to  compute  the  K's  and  solve  for  these  parameters  again.  This 
gives  the  second  approximation.  Repeat  the  cycle  until  the  desired 
accuracy  is  obtained.  To  make  the  convergence  more  rapid,  after  the 
third  approximation,  instead  of  following  the  above -outlined  method 
rigidly,  one  may  extrapolate  the  results  and  use  the  values  of  the  A's 
so  determined  to  repeat  the  process.  This  alternative  method  is  found 
to  be  especially  time  saving  in  the  computation  for  the  cases  where  the 
Mach  number  is  relatively  high.  A  detailed  example  of  the  computation 
may  be  found  in  reference  4. 


The  parameters  A±j/a0  for  flows  at  various  undisturbed- stream 

Mach  numbers  are  computed  and  are  tabulated  in  table  I.  In  order  to 

see  the  convergence  of  the  series,  these  parameters  are  computed  by 

taking  one  parameter  only,  two  parameters,  three  parameters,  and  so 

forth,  up  to  six  parameters.  The  maximum  velocities  at  various  Mach 

numbers  are  tabulated  in  table  II  for  all  the  six  cases.  The  convergence 

of  the  series  is  seen  to  be  satisfactory  except  at  Mq  =  0.5  where  the 

series  is  apparently  divergent,  and  with  four  terms  the  maximum  velocity 

is  already  greater  than  the  maximum  allowable  velocity  a  .  For  y  -  2 

iii  ax 

(% iax\2  2 

— j  =  2  +  MQ  .  The  velocity  distribution  over  the  cylinder  at 


Mq  =  0.4  computed  by  the  variational  method  is  plotted  in  figure  2  as 
are  the  results  by  the  Rayleigh- Janzen  method  to  the  third  approximation 
as  computed  by  Imai  (see  reference  6).  At  various  Mach  numbers  of  the 
undisturbed  flow,  the  maximum  local  Mach  number  is  plotted  in  figure  3 
and  the  local  Mach  number  over  the  surface  of  the  cylinder  is  plotted 
in  figure  4.  The  velocities  computed  by  the  Rayleigh- Janzen  and  the 
present  method  show  excellent  agreement. 


In  order  to  study  the  effect  on  the  velocity  of  using  7=2,  the 
next  value  closer  to  the  isentropic  one,  that  is,  7  =1.5,  is  now  taken 
to  compute  the  velocity.  Taking  only  the  first  term  in  the  velocity- 
potential  series,  the  equation  bl/bA^  =  0  is  obtained  as  follows: 

(2-1Ull'W2u3  -  3.28889U5)  -  * 

5.92804Ulf)A11  -  (3.21005U2  -  0.82222qfflax2)lJA112  -  (l.0^66v2  - 
0.10159qmax2)A113  -  0.l6490UA114  -  0.010564An5  -  o  (28) 


Ik 
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The  maximum  velocity  with  7  =  1.5  is  also  included  in  table  II  for 
comparison.  It  is  seen  that  these  velocities  are  smaller  than  the  cor¬ 
responding  ones  computed  by  taking  7=2.  However,  the  velocities  also 
increase  with  increased  number  of  terms  in  the  f  series.  Since  only  a 
finite  number  of  terms  can  be  taken  in  the  numerical  calculation,  the 
velocities  computed  will  always  be  smaller  than  the  corresponding  exact 
values.  In  such  cases,  the  effect  of  taking  7=2  is  to  increase  the 
velocity,  which  compensates  to  some  extent  the  effect  of  taking  only  a 
finite  number  of  terms. 

The  pressure  coefficient  at  a  point  whose  pressure  is  p  may  be 
defined  as  follows: 


The  pressure  coefficient  is  computed  from  the  above  formula  and  is 
plotted  in  figure  5.  Results  by  the  Rayleigh- Ja nzen  method  and  the 
variational  method  compare  very  well,  with  the  former  yielding  slightly 
lower  results  at  small  angles.  The  pressure  coefficient  by  the  Rayleigh- 
Janzen  method  was  computed  by  taking  7  =  1.405  in  equation  (29),  which 
becomes  upon  substitution: 


Both  the  variational  and  Rayleigh- Janzen  pressure-coefficient  values  are 
lower  than  those  obtained  by  the  Von  Karmdfn-Tsien  method  and  generally 
higher  than  the  Prandtl-Glauert  approximation.  For  further  comparison, 
the  pres sure -coefficient  ratio  is  plotted  in  figure  6  where 

Cp^  is  the  pressure  coefficient  of  the  corresponding  incompressible 
flow. 
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FLOW  PAST  A  CURVED  SURFACE 


In  order  to  obtain  further  evidence  of  the  applicability  of  the 
method,  the  compressible  flow  past  a  thin  curved  surface  will  next  be 
considered.  This  example  is  carried  out  for  two  purposes.  First,  from 
the  previous  example,  the  variational  method  is  seen  to  give  very  good 
results  for  flow  past  a  thick  body  where  the  Mach  number  is  relatively 
lowj  it  is  to  be  ascertained  now  whether  the  method  will  also  give  good 
results  for  flow  past  a  thin  body  at  high  Mach  numbers.  Second,  there 
may  be  some  difficulty  in  carrying  out  the  method  when  an  arbitrary 
mapping  function  is  used.  If  this  is  so,  a  study  should  be  made  to 
remove  such  difficulty. 

The  compressible  flow  past  a  thin  bump  has  been  carried  out  by 
Kaplan  (reference  1)  by  the  perturbation  method  in  terms  of  a  thickness 
parameter.  Results  are  given  for  flow  past  a  series  of  bumps  of  various 
thicknesses.  Kaplan's  bump  may  be  obtained  by  mapping  conformally  a 
circle  of  unit  radius  in  the  complex  z-plane  into  the  complex  {[ -plane  by 
the  mapping  function 


,  ^2  ,2 
5  =  z  + 

2  3z3 


(32) 


Again  take  six  parameters  in  (f  as  in  the  previous  example,  noting 
that  z  «  re*^.  The  equations  =  0  may  now  be  written  in  the 

following  form: 

*N  ->v  \  ~  2lC  >. 


max 


Mil 


ieSAjJ.st  d9  T  ~ 


r  Cn) 


de 


r  dr 


p  n  2jt 
^1  ^0 


M  a  M  M 

j 


dr 

r 


obtained  from  the  mapping  function  (32)-  The  quantity  6^  -  1  when 
i  =  j  =  1  and  5ij  =  0  for  other  values  of  i  and  j. 

Substituting  into  the  integrals  (33)  the  expressions  (25)  and 
carrying  out  the  integration  first  with  respect  to  the  angle  0,  the 
resulting  integrals  are  of  the  type 

,  •  n  2lt 

In  =  /  Xf(sin  0  cos  0)d0 
J  0 

This  integration  can  be  somewhat  simplified  by  the  use  of  the  theory  of 

residue.  Let  Z  =  ei0  and  express  the  integrand  of  I  as  a  function 
of  Z.  Each  I  then  consists  of  a  sum  of  integrals  of  the  two 
following  types:  -  ■ 
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There  are  altogether  66  In  integrals  encountered  in  the  calcula¬ 
tion,  and  they  are  given  in  appendix  A.  With  these  expressions  of  In 
substituted  into  the  integrals,  the  resulting  integrands  are  rational 
functions  of  r  with  coefficients  in  d.  It  is  observed  that  the  final 
coefficients  of  the  unknown  parameters  A^j  are  sums  of  the  following 

functions: 


(rh  *  d2)^  -  ^(r*  -  l) 


where  n  may  be  either  positive  or  negative.  There  are  altogether  20 
such  Yn  expressions  encountered  in  the  calculation,  and  they  are 

given  in  appendix  B.  It  may  be  noted  that  terms  like  [loge  (1  ~  r^7j  r 

and  (loge  r )r_00  appear  after  the  integration.  These  terms  however 

will  be  canceled  by  each  other  in  the  final  summation.  This  fact  may 
conveniently  be  regarded  as  a  check  to  the  correctness  of  the  computati 


The  numerical  values  of  d2  may  be  substituted  at  this  stage  befo: 
the  final  integration  with  respect  to  r  is  carried  out;  otherwise  the 
computation  will  be  extremely  laborious.  The  integrands  now  may  be  tak 
as  the  products  of  r  expressions  and  In  and  then  the  final  answers 
can  be  immediately  written  down  from  the  known  Yn  values.  There  are 
eight  r  expressions  involved,  and  they  are 


P 


R 


K 

M 


3r 

1  - 


_1 

r^ 


1 

r 

Q-i-1, 

1 

rH 

i— 1 

1 

Cl 

■v 

5  7 

5r  7r f 

i 

r  _  __L  +  JL 

r2 

L  2  +  4 

1 

As  a  numerical  example,  d2  =  0.075  is  taken,  which  corresponds  to 
a  bump  with  a  maximum  thickness  of  5*13  percent  of  the  chord.  The  shape 
of  the  bump  is  shown  in  figure  7-  For  d  =  0.075>  the  values  of  Yn 
are  calculated  and  included  in  appendix  B.  After  dropping  a  common 
factor  of  4,  integrals  (33),  which  are  equivalent  to 


cd/aAil  =  0  Si/6A13  =0  Sl/dA31  =  0 

6i/SA33  =  o  6i/6a15  =  o  6i/6A51  =  o 
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"become 


d.5185190^2  -  i.193255U2)a11  +  o.55j+586u2A13  +  ^0. 0814815^ 
.i69149U2)a31  +  0.0285596u2a33  +  o.216663U2a15  +  ^0.0285714^. 
.0546889U2)a51  =  % 


0.554586u2An  +  (3- 333333^^  -  8.188034u2)a13  +  0.112085U2A31  + 
^0.400000qmax2  -  i.012468u2)a33  +  2 . 347577U2a15  +  0.0457158u2a51  =  Kg 

(o.08l48l5qfflax2  -  0.169149U2)a1;l  +  0.112085U2A13  +  (o.0225l85qmax2  - 
0.0420399U2)a31  +  0.00795193U2A33  +  0.0706980U2A15  + 

(o.0102858qmax2  -  0.0178331U2)a51  =  K3 

o.0285596u2Au  +  (o.4oooooqmax2  -  i.012468u2]ai3  +  o.oo795193U2a3i  + 
(o.o693333gmax2  -  0.i676i6u2)a33  +  0.3198i8u2a15  +  0.00397346u2a51  = 


0.216663u2ai;l  +  2.347577U2A13  +  0 . 0706980U2A31  +  0 . 3l98l8u2A33  + 
(8.962963qmax2  -  22.05199U2)a15  +  0.0356363U2a51  =  K5 

^0.02897i4qrripy2  -  0.0546889U2 |ali  +  o.o457158u2a13  + 

(0. 0i02857qmax2  -  ,0.0178331U2)a31  +  0.00397346u2a33  + 


0.0356363U2A15  +  (0.0054577qmax2  -  0.oo887571U2)a51  =  Kg 


The  Kn  expressions  are  found  in  appendix  C. 
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These  equations  (3^)  are  solved  again  hy  the  method  of  successive 
approximation  as  outlined  in  the  previous  section.  At  relatively  high 
Mach  numbers,  however,  the  method  as  outlined  before  becomes  oscilla- 
tively  divergent.  In  such  cases  the  average  values  of  A±^  between  two 
consecutive  cycles  of  computations  may  be  used  for  the  next  approxima¬ 
tion.  The  process  is  then  again  convergent.  Such  a  method  was  also 
found  necessary  by  Wang  in  a  different  problem  where  another  set  of 
cubic  simultaneous  equations  was  involved  (reference  7)* 

At  Mq  =  0.9,  the  method  described  above  again  failed  when  six 
parameters  were  used.  Many  schemes  of  computation  were  tried  and  it 
was  not  possible  to  calculate  the  values  of  the  parameters  by  methods 
of  successive  approximation.  It  was  possible,  however,  to  calculate 
these  values  when  up  to  four  parameters  were  taken.  This  indicates  that 
at  this  Mach  number  irrotational  solution  may  not  exist  (reference  8) 
and  shock  waves  have  possibly  occurred.  It  can  also  be  seen  that  the 
velocity  computed  by  Kaplan  begins  to  diverge  at  about  Mo  =  0.9.  The 
velocity  series  computed  by  Kaplan  becomes  in  this  case 


Q  2  3 

=  1  +  a.jt  +  a^t  +  a^t°  +  •  •  • 

*  1  +  0.176472  +  0.050675  +  0.059108  +  .  .  . 

where  t  Is  the  maximum  thickness  ratio  and  is  equal  to  0.051282  in 
this  case. 

The  values  of  k±^a.Q  are  given  in  table  III  at  Mq  =  0..5>  0.75 > 

0.83>  and  0.90.  The  velocity  distribution  over  the  surface  of  the 
bump  is  given  in  table  IV  and  plotted  in  figure  7  to  compare  that  found 
by  the  variational  method  with  the  results  by  Kaplan.  The  agreement  is 
very  good.  At  Mq  =  0.9>  the  maximum  velocities  q/U  with  one,  two, 

three,  and  four  parameters  are  1.1179*  1.2628,  1.2716,  and  I.2636, 

respectively.  The  successive  values  of  the  velocity  do  not  show  a  sign 
of  divergence  in  this  case.  The  velocities  at  various  points  on  the 
surface  of  the  bump  are  obtained  as  follows.  For  the  mapping  func¬ 
tion  (32),  the  shape  of  the  bump  in  the  £  -plane  is  given  by  the  following 
parametric  equations. 


1=2  cos  G  -  “~(3  cos  G  -  cos  3$) 
T)  =  ~-(3  sin  G  -  sin  30 ) 


N 

* 

) 


(35) 
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The  length  of  chord  is  given  by  (26)0=Q  or  ^2  -  and  the  maxi¬ 
mum  thickness  is  given  by  (2r| )  .  or  4d2/3.  The  thickness  ratio  t 

p  ^  v7 — Jt  j  2 

is  then  equal  to  2d  /(3  -  d2).  If  the  unit  length  is  taken  as  the 
semi  chord  of  the  shape,  equation  (35)  "becomes 

t,  ^ 

X  =  cos  6  -  jj-(cos  6  -  cos  3^) 

r  {36) 

Y  =  ^(3  sin  0  -  sin  30)  J 

For  different  values  of  X,  9  may  be  found  from  equations  (36).  On 
the  surface  of  the  bump, 

-  (1  -  d2)cos  20  -  d2  cos  +  £(1  -  d2)sin  20  + 

d2  sin  (37) 


The  velocity  is  then  given  by  the  following  expression: 


,r  Se/r=l 

“ s 

dz 


1  [7  0 . 66667 A11  0.133333A31  0.05711*2A51\ 

[f  aoMo  +  +  7 


sin  0  + 


/2A13  o.4a33\  3.33333A15 

W +  ^%rln  39  +  “^7^ sin  59J  (38) 

For  different  values  of  X,  the  corresponding  values  of  sin  G ,  sin  3 Q  y 

d.  C 

sin  50,  and  —  are  given  in  table  V. 

The  pressure  coefficients  at  various  Mach  numbers  are  given  in 
table  VI  and  are  plotted  in  figure  8.  The  pressure  coefficients  com¬ 
puted  by  the  variational,  the  perturbation,  the  Von  Kahnffia-Tsien,  and 
the  Prandtl-Glauert  methods  are  given  in  table  VII  and  are  plotted  in 
figure  9  for  comparison. 
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RESULTS  AND  DISCUSSION 


Compressible  flows  past  a  circular  cylinder  and  a  thin  curved  sur¬ 
face  have  been  computed  by  the  variational  method.  In  the  case  of  flow 
past  a  circular  cylinder,  the  velocities  and  pressure  coefficients  com¬ 
puted  are  compared  with  the  results  obtained  by  Imai  who  used  the 
Rayleigh-Janzen  method  up  to  the  third  approximation.  In  the  case  of 
flow  past  a  thin  curved  surface,  the  velocities  and  pressure  coeffi¬ 
cients  are  compared  with  the  results  obtained  by  Kaplan  who  used  the 
perturbation  method  in  terms  of  a  thickness  parameter  up  to  the  third 
approximation.  The  agreement  of  the  results  between  the  present  method 
and  the  existing  methods  is  very  good.  It  thus  shows  that  the  varia¬ 
tional  method  will  give  good  approximate  solution  both  for  thick  and 
thin  bodies  and  both  at  low  and  high  Mach  numbers. 

In  carrying  out  the  variational  method,  7  is  taken  to  be  2  instead 
of  the  isentropic  value  which  for  air  is  1.405.  A  study  is  made  by 
taking  7=1.5  in  the  case  of  flow  past  a  circular  cylinder.  The 
results  show  that  the  effect  of  using  a  larger  value  of  7  is  to 
increase  the  maximum  velocity  slightly.  The  same  conclusion  may  be 
obtained  by  a  study  of  the  results  of  the  Rayleigh-Janzen  method. 
Theoretically,  if  one  can  take  a  complete  set  of  functions  in  the 
velocity-potential  series,  exact  solution  may  be  obtained.  In  actually' 
carrying  out  the  method,  however,  only  a  finite  number  of  terms  can  be 
taken.  Although  the  introduction  of  more  terms  into  the  solution  does 
not  involve  any  additional  mathematical  difficulties,  each  new  term 
appreciably  increases  the  computational  labor.  With  six  terms,  it  is 
observed  that  the  maximum  velocities  increase  with  an  increased  number 
of  terms.  Thus  the  slightly  higher  values  for  maximum  velocities  are 
increments  in  the  right  direction.  By  using  a  modified  pres sure -density 
relationship  instead  of  the  isentropic  one,  the  pressure  coefficients 
are  hardly  affected  by  an  increase  in  7. 

When  the  undisturbed  velocity  increases  to  a  value  at  which  shock 
waves  will  probably  occur,  the  nonlinear  terms  of  the  undetermined 
parameters  are  no  longer  small  compared  with  the  linear  terms.  The 
method  of  successive  approximation  in  solving  these  simultaneous  equa¬ 
tions  suddenly  fails.  None  of  the  many  schemes  tried  worked  when  more 
than  four  parameters  were  taken.  An  investigation  has  been  carried 
out  to  study  the  reason  why  it  is  so.  It  was  found  (reference  8)  that 
the  irrotational  flow  solution  fails  to  exist  when  the  Mach  number  is 
increased  beyond  a  certain  limiting  value.  For  the  cases  computed, 
it  appears  that  this  limiting  value  decreases  as  the  number  of  terms 
increases.  This  explains  why  the  method  of  successive  approximation 
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fails  to  work  with  more  parameters  while  it  was  still  possible  to 
obtain  a  solution  with  fewer  parameters. 
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I^=  X  COS2@  COS23^  d0 


4(1  +  d2)(r4  +  d2)(r4  -  d4)  (r4 


[7(1  +  d2)r12  +  2(1  -  d4)r10  + 


(1  +  d*)(l  -  3d2  +  d4)r8  +  2(1  -  d°)rb  +  (l  -  d4)(l  -  db): 


2d4(l  -  d4)r2  -  d4(l  +  d6 


Ic  =  /  X  cos  0  cos  3 Q  cos  50  d0 
J  0 


1  +  d2)rl4  +  2(1  -  d4)r12  + 


4(l  +  d2)!1-^  +  d2)(r4  -  d4)(r4-  l)r2  LA-  ^ 

(l  +  d2)(l  -  d2)2r10  +  (l  -  d8)r8  +  (2  -  d4  -  d6  +  2d10)r6  + 
(l  -  d4)  (l  +  d8)r4  -  2d4(l  +  d8)r2  -  d4(l  -  d8) 


-  /  X  COS20  COS25^  <10 


4(1  +  d2)(r4  +  d2)(r4  -  d^)  (r4  - 


2(1  +  d2)r16  +  2(1  - 


d4  r14 


2d2^l  +  d2)r12  +  (l  +  d10)r8  +  2^1  -  d12)r 


l4)(l  -  d10)  r4  -  2d4(l  -  d8)  r2  -  d4(l  +  d10 


4 
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n2ic 

*10  =  ^  cos  ^  cos^50  d@ 

H1  +  d2)(rl*  *  4*  -  ^  ^ + 4  ^  ' 4 

3d^(l  +  d2)r12  -  3d4(l  -  d4)r10  +  (l  -  dl6)r6  +  (l  +  dl8)r4  - 
d4(l  -  d12)r2  -  d4(l  +  dl4)j 


^11  =  I  ^  COS20  d0 


jtr 


^(l  +  d2)(r4  +  d2)(r4  -  d^)(r^  -  l) 

(l  +  d2)(l  +  d^)^  +  d^(l+dfl 


(l  +  d2)r8  - 


-12 


2n 


\  si el20  cos23^  d0 


4(l  +  d2)(r^  +  d2)(^-d^)(r^-.l)L"(- 
(l  +  d2)(l  +  d2  +  d4)r8  +  2(1  -  d8)r6  - 
2d4(l  -  d4)r2  +  d4(l  +  d6) 


2(1  +  d2)r12  -  2^1  -  d^)r10 


1  -  a4)(l  -  d6)r4  - 


14 
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1-^3  =  J  X  sin20  cos25 0  d 6 


2d2^l  +  d2)r^2  -  (l  +  d'*'^)r®  +  2^1  -  d^2)r8  -  ( 

0,W,  ^  2  ,lo\] 

2d  11  -  d  }v  +  d  ^1  +  d  )\ 


1  +  d2)r16  -  2(1  -  dh)rlk 


1  -  d^)(l  -  d10^ 


=  J  X  sin  0  cos  6  cos  30  d0 


4(l  +  d2)(r^  +  d2j(r^  -  dk)(rk  -  l) 


1  +  d2  r10  -  1  -  dMr8 


(l  +  d2)(l  +  d4)r6  +  (l  -  d8)vk  +  d4(l  +  d2)r2  -  d4(l  -  dk 


-j  X  sin  6  cos  6  cos  5 0  d0 


4(l  +  d2)(ri|’  +  d2jfr^  -  dMlr^  - 


1  -  dl)r10  -  (1  +  a6)r8 


(l  -  d8)rS  *  (l  +  d^)(l  +  d^r1*  *  d4(l  -  d^r2  -  a*(l  ♦  d«)_ 


28 
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1  +  d10  r4  +  d4  1  +  d 


32 
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p  2jt 

1^2  =  J  X  sin2©  sin25S  d© 


38 
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n  2rt 

1^9  =  J  X  sin2  3 0  cos  30  cos  5 6  d 0 


Ip-,  =  /  X  sin  3@  sin  5^  cos  3^  cos  50  <3-0 

J0  . 
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X  cosJ3$  cos  50  d 6 


_  7T _ 


+ 


(l  +  a^r3-1*'  +  (3  -  -  d6  +  3d10)r10  -  3d^(l  +  a6)r6  + 

(1  -  al6)r^  -  a^(i  -  a12) 
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APPENDIX  B 


THE  Yn  EXPRESSIONS 

After  integration  the  functions  Yn  can  be  put  into  closed  forms 
in  terms  of  tan'1  d’1  and  loge  ^i!,  vhere  lofe  lndlcate,  natural 

eS^-J^eT^dilSenS1”;  °f,  d'  the  'Station 

nite^series  £ 

follows:  d-  In  SUCh  f0rmS^  the  Unctions  Yn  are  as 


*1  = 


1 _  fD  1 

2(1  +  d2)2(l  -  d2H2  ”  2  loge  2  +  1  -333333333  -  0 . 200000000a2  + 


0. 3428571 W-  -  0.053968253d6  +  0. 202020202d8  -  0 . 024975025a10 
0.143589743d12  -  0.0l4379085d11^  +  0.111455107a16  +  .  .  \ 


Y3  =  - 


4  (i  +  d2)2(l  -  d2)^1  +  l0Se  2  +  1-  000000000  -  0  •  5°0°00000d2 
'•833333333d  -  0.083333333d6  +  0.450000000a8  -  o.033333333a1( 
•3°9523809d12  -  0.017857142a14  +  0.236111111a16  +  .  \ 


Y5  =  ’  2(1  +  a2)  2 (1  .  a2)  ff  ‘  ^  l0Se  2  +  ^O00000000  -  0.333333333d2  + 


0.533333333d  -  0.009523809d6  +  0.253968253d8  -  0 . 002020202d10  + 
0.167832167a12  -  0.000732601a14  +  0.125490196a16  +  .  .  .) 
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(-D  -  log, 


2  +  dr  -  1 


.  500000000d1)'  -  0.1666666674® 


“(l  *  i2)2(l  -  d2)V 

0.583333333d®  -  0.050000000d10  -  o.366666667d12  -.0.0238095254 
0. 2678571 ^Sd1® 


2(1  +  d2)2  t-  -  d2) 


L»  +  1 

^2  2 


log  2  +  d2  -  1.333333333d 


6  .  0 . 3^28571 42d®  -  0.079365080d10  -  0.2020202024' 


0.466666667a 


032167833d114'  -  0. 143589743d1®  +  . 


l(l  +  d2)2  (1  -  d2) 


D  +  loge  2  +  a 


^  -  1.5000000004® 


8  +  0.4l6666667d10  +  0.4500000004' 


0.833333333d 


0.133333333d14  +  0 . 309523809dl6  + 


_  _ 1 _  .  (p  +  log  2  +  1.500000000  -  0.333333333d  + 

4(l  +  d2) 2  (l  -  d2)  6  ^ 

0.583333333d4  -  0.ll6666667d6  +  0.3666666674®  -  0.0595238094  + 

0.2'67857l42d12  -  0.036lllllldl4  +  0 . 2111111Hdl6  +  .  .  .) 


1 _ {H  .  1  log  2  +  1.533333333  -  0.l42857l42d  + 

/  p\  2  /  \2  2  e 

2(l  +  d2)  (l  -  d2) ' 

0.253968253d14  -  0. 0591630604®  +  0.1678321674®  -  0.032478632410  + 
0. 125490 196d12  -  0 . 0205460l6dl4  +  0.100250625d  +  •  •  •) 
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Y-5  =-77 - IT§7 - rr(o  +  loge  2  +  1.833333333  -  0.250000000d2  + 

b[i  +  a2}  [1  -  d2) v 

0.4500000001 4  -  0 . Il6666667d^  +  0 . 309523809d8  -  0 . 067857142110  + 
0.236lllllld12  -  0.044444444dl4  +  0.19090909H18  +  .  .  .j 

y_7  =  -  -7^ - ^y(|  -  |  ioge  2  +  1.676190^75  -  o.iiiiiinid2  + 

0.20202020214  -  0 . 056721056d6  +  0 . 143589743d8  -  0.03458ll05d10  + 
0.11l455107d12  -  0.023327549114  +  0.091097307dl8  +  .  . 

y-9  - - 7 - 77^7 - r(D  +  loge  2  +  2.083333333  -  0.200000000d2  + 

4(1  +  d2]  (1  -  d2)v 

0.366666667d4  -  0.10952380918  +  0.267857l42d8  -  o.o69444444d10  + 
0.2liiiiuid12  -  0.04805l949d14  +  o . I74242424dl6  +  .  .  .) 

Y-n  =  — — - 7757 - rr(!  -  7  log  2  +  1.787301587  -  0.0909090911 2  + 

2(1  +  d2)  (1  -  d2)' 

0.l67832l67d4  -  0 . 05268o652d8  +  0.125490196d8  -  0.034532031110  + 
0.100250625d12  -  0.024430640dl4  +  0.083478260d18  +  .  .  .j 

y-13  =  -  — - —77 - t(d  +  log  2  +  2.283333333  -  0.i66666667d2  + 

4(i  +  d2)  (1  -  d2)  e 

0.309523809d4  -  0.101190475d6  +  0.236lllllld8  -  0.068253969d10  + 
0.l9090909ld12  -  0. 04924242 4dl4  +  0 . l6o256409dl6  +  .  .  .) 
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-15  = 


"  “I - T51 - Ttl  “  I  loSP  2  +  1-878210678  -  0.076923076d^  + 

2(1  +  a2)2(i  -  a2)^2 

0.143589743a1''  -  0.048567120a6  +  0.111455107a6  -  0.033583958a10  + 


0.091097307a 


12  -  0.024654731a11  +  q. 077037037a16  + 


1  +  a‘ 


)2(i  -  *2) 


D  +  log.  2  +  2.450000000  -  0.142857142a  + 


0.267857142a1  -  0.093253969a6  +  0.211111111a8  -  0.065909091a10  + 
0.174242424a12  -  0.049145298a11''  +  0.l4835l647al6  +  .  .  .) 


- - - T5-7 - r(§  -  §  i°Se  2  +  1-955133754-  -  0.066666667a  + 

2(1  *  d2)2(i  -  2 

0.125490196a1  -  0. o4478844oa6  +  0.100250625a8  -  0.0 32273778al°  + 


0.083478260a12  -  0.024405459a11  +  0.071519795a16  +...), 


4(1  +  a2)2(i  -  a2) 


d  +  logQ  2  +  2.592857142  -  0.125000000a  + 


0.236111111a1  -  0.086111111a6  +  0.190909091a8  -  0.063131313a10  + 
0.160256409a12  -  0.048351647a11  +  0.138095238a16  +  .  .  .) 


Y_23  =  - 


(1  +  a2f(i  -  a2) 


|  -  |  loge  2  +  2.021800421  -  0.058823529a2  + 


0.111455107a1  -  0.041427096a6  +  0.091097307a0  -  0.030846682a  + 

0.077037037a12  -  0.023900748a11  +  0.066740822a16  +  .  . 
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Y  2  =  -  — - -|t - t(d  +  loge  2  +  2.717857142  -  O.llllllllld2  + 

4(l  +  d2J  (l  -  d2]V 

0 . 2lllllllid^  -  0.07979T980d6  +  0.174242424d8  -  0.060256409d10  + 
0.l4835l647d12  -  0.047l86l47dllt'  +  0.129166667a16  +  .  .  .) 

Y  27  =  -  -7 - — r^7 - -r(i  -  \  log  2  +  2.080623951  -  0.052631578d2  + 

1  2(l  +  d2r(l  -  d2)^ 

0.100250625d1)-  -  0.038465729d6  +  0.083478260d8  -  0.0294l7990d10' + 
0.07i5i9795d12  -  0.023262562a111  +  0.06256i094dl6  +  .  .  .J 

where 


E  =  (loge)r=oo 

For  d2  =  0.075>  the  values  of  Yn  are  as  follows: 

Y±1  =  -0.163282840  -  D’  +  E 

Y9  =  0.193584561  -  D’ 

Y7  =  -0.146562469  -  D’ 

Y5  =  -0.295349926  -  D’ 

Y3  =  -0.388303682  -  D’ 

Yx  =  -0.455432461  -  D’ 

Y_!  =  -O.5O783H36  -  D' 

Y_3  =  -0.550751149  -  D' 

Y_5  =  -0.587074905  -  D' 

Y_t  =  -0.618549324  -  D' 


y_9  =  -0.646311118  -  D’ 

Y_u  =  -0. 671140539  -  D' 

Y_13  =  -0.693595951  -  D' 

Y  -,c  =  -0.714090490  -  D' 
-15 

Y_i7  =  -0.732938314  -  D’ 
Y_i9  =  -O.75038382I  -  D' 
Y-21  =  "°* T66620856  -  D' 

Y_23  =  -O.78I805725  -  D' 

Y_25  =  -0.796066296  -  D' 
Y_2T  =  -0.809508489  -  D' 


0 


.  233873678 [loge( 


1  -  r‘ 
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APPENDIX  C 

THE  Kjj  EXPRESSIONS 

The  Kn  expressions  encountered  in  the  calculation  of  flow  past  a 
curved  surface  by  the  variational  method  are  as  follows: 

%  =  0.059999539U3  +  0.i8295it90UA112  -  0.2890247071%-%  3  -  o.099215517UA11a15  + 

0.0395^924toA11A31  -  0. 0551459 4ouai;la3 3  +  o.oii687o44ua11a51  +  1. 52197  Vk>7UA132  - 
1.708453568oa13a15  -  o.039032242UA13a31  +  0.4623i8i73UA13A33  -  o.oio547763UA13a51  + 

4 . 0627969761% j2  -  o.ooo252462UA15a31  -  0.203909l65UA15A33  +  0.0009786001%^%  - 
o.ooo3874o6ua312  -  0.015012518ua31a33  +  o.oo380236oua31a51  +  0.037053721UA332  - 
0 . oo8i34i89UA33%1  -  0. 0007491111%  j2  +  0.033220706ai;l3  +  o.022268910Ai;l2a13  + 
0.0l6688092A112AJ5  +  0.023629862A112A31  -  o:oo4191746a112A33  +  0.010268355A112A51  + 
0.586559872AnA132  -  0.052992785A11A13A15  +  0.030572932A11A13A31  + 
o.196566992A11a13a33  +  0 .0178i6638a11a13A51  +  x.590352o68a11a152  + 
0.027583320A11A15A31  +  o.026348o6xa11a15a33  +  o.Oi6797344a11A15A51  + 
0.007609784a11a312  +  o.oo8ii84o6a11a31a33  +  0.007721948a11a31a51  + 
o.oi94o493XA11A332  +  0.00503720XA11a33a51  +  0 . 002i47644a11A512  +  o.35534846oa133  + 
1.998793077A132A15  +  o.126029525A132a31  +  0.192711803A132A33  +  o.0545548o4a132a51  + 
2.546743570A13a152  +  0.093686l39A13A15A31  +  0.725826526a13A15A33  + 
o.o66i65622a13a15a51  +  o,oo7521368a13a312  +  o.o496o6oi4a13a31a33  + 
o.oo8865o88a13a31a51  +  o.038026476a13a332  +  0.021958499A13a33a51  + 

0.002747290A13A512  +  1.182535859A153  +  0 . 354628l07A152A31  +  0.493373000A152A33  + 
0.1588845i8a152a51  +  0.006593000A15A312  +  o.02623556oa15a31a33  + 
o.oo8o96717A15a31a51  +  o.o68307i46a15a332  +  0.017553086a15a33a51  + 

0.002597463A15a5x2  +  o.ooio63o6oa313  +  0.002385806a312a33  +  O.ooiSiSgiOA^A^  + 
0.00549159XA31a332  +  0.003005606a31a33a51  +  o.ooiio698oa31a512  +  0.002670255A333  + 
0.0028i6774a332a51  +  o.ooo98738oa33a512  +  o.ooo237098a513 


=  -  0.130908999U3  -  o.i44512351»ua112  +  3.o4391t88ltoA11A13  -  1 . T08453568UA11A15  - 

0.039032242UA11a31  +  0.4623i8i73UAXiA33  -  0.010547763UA11a51  +  3.21377h8ua132  + 
ii.857940352UA13a15  +  0.346476o64ua13a31  +  i.031279To4ua13a33  + 
o.091769946ua13A51  +  7-73l27342lUA152  -  0.3498o8o85UA15A31  +  1.726o43683UA15a33  - 
0 .  i46o59989UA1jA^1  +  o.ooo692852UA312  +  o.058oo4535ua3^a33  +  0.002368c>92UA3]_A^  + 
0.087200038ua332  +  0.0609306001^3^!  +  0.001909217UA5-,2  +  0.007422970A1;L3  + 
0.586559872Ai;l2A13  -  0.026496392AU2A15  +  0.015286466ai;l2A31  + 
o.098283496a112a33  +  0.008908319A112A51  +  1.o66o45380A11A132  + 
3-997586i54a11a13a15  +  0.252059050A11a13a31  +  o.3854236o6a11a13a33  + 

0 . 109109680A11A13A51  +  2.564743570A11a152  +  0 . 093686l39Al;LA15A31  + 

0.725826526a11a15A33  +  0.06616^622^^^  +  0.007521368ai;la312  + 
0,o496o6oi4a11A31A33  +  o.oo8865088A11A31A51  +  0.038o26476aua332  + 
0.021958499Ai;lA33A51  +  0.002747290A11A512  +  3.609289260^  +  8.71963914oa132a15  + 
0.208248666A132A31  +  ij«i1741922A132A33  +  0.085998486a132A51  + 

25.4l42i8o49A13A152  +  0.920228170A13A15A31  -  3.000980028a13a15A33  + 

o.4551o4944a13a15a51  +  0.039746i67A13A312  +  0. 086959 326a13a31a33  + 

o.o42256i62A13a31a51  +  o.324o78519A13a332  +  o.o4ox62i46a13a33a51  + 
0.012506092A13A512  +  10.871460415A153  +  o.504585238a152a31  +  4.344o47335A152a33  + 
0.210736557A152A51  +  o.034424058a15a312  +  0.i68542ll6A15A31A33  + 

0.04555l46lA15A31A51  +  0 . 28989 4836a1^a332  +  0.040344011A15A33A51  + 

o.01564565IA15a512  +  0.001033455A313  +  o.Oo8o85122A312A33  +  o.ooi92735lA312A51  + 

0.009687255a31a332  +  o.oo9037828a31a33a51  +  o.oou82158a31a512  +  0.020369802A333  + 

I 

o.oo4917098a332a51  +  o.oo2764734a33a512  +  o.ooo247234a513 
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K3  =  0.009387226U3  +  0.019774622UA1;L2  -  0.0390322teUA11A13  -  0 . 000252462UAuA15  - 

o.ooo7748i2ua1:la31  -  0 . 015012518uai;la33  +  o.oo38o236oual1a^1  +  o.173238o32ua132  - 
o.3498o8o85UA13A15  +  0.001385705UA13A31  +  o.058oo4535UA13A33  +  o.od2368o92UA13A51  + 
oA27873864ua152  +  o.020248502ua15a31  -  0. 03594720 3UA15a33  +  0.007519578ua15a51  - 
o.oo253184oua312  -  0 . oo4444278ua31a33  -  0.000566987UA31a51  +  0.003749007UA332  - 
0.0020907^90A33A51  -  0.00107891hua512  +  o.O0787662OAn3  +  0.015286466a112a13  + 
0.013791660A112A15  +  0.00760978Uai;l2a31  +  o.oo4o59203Ai:l2a33  +  0.00386097^^  + 
0.126029525A11A132  +  0.093686139A11A13A15  +  o.oi5042736a11a13a31  + 

0 . 04960601  i+a11  A-j^  3a33  +  o.oo8865088a11a13a51  +  0.35 4628io7a11  a^2  + 
o.oi3186oooai;la15a31  +  o.o2623556oa11a15a33  +  o.oo8o96717A11a15a51  + 
0.00318918oai;la312  +  0.004771612^^^  +  0.003633220Ai;la31a51  + 

0.005491591 a11a332  +  0.003O056o6a11a33a^1  +  o.ooho698oa1;1a512  +  o.o694i6222A133  + 

0.46oii4o85A132a15  +  0.039746i67A132a31  +  o.o4347-9663A132a33  +  0.02112808^^  + 
0.504585238ai3A152  +  0.0688481^^^  +  0.1685421^^^  + 
o.o4555i46ia13a15a51  +  0:003100365A13a312  +  0.0i6i70245A13A31A33  + 
0.003854703A13a31a51  +  0.009687255A13a332  +  0.0O9037828a13a33a51  + 
o.ooii82158a13a512  +  o.23o447196a153  +  o.12o6oo559a152a31  -  o.26o365077A152a33  + 
o.o67340386a152a51  +  0. 00270887^^^  +  0.01538i835a15a31a33  + 
o.oo329776ia15a31a51  +  o.oi764i290A15A332  +  o.oioo744i5a15a33a51  +  ' 
0.001024491A1^12  +  0.00055l450A313  +  0.0013234l0A312A33  +  0.001028092A312A51  + 
0.002235226a31a332  +  o.ooi7o646ia31a33a51  +  o.ooo672957A31a512  +  o.ooo845763A333  + 

0.00134l68lA332A51  +  0.0005698l5A33A512  +  0.000153120A^13 
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=  -0.016023195U3  -  o.027572970UA112  +  0 . 4623i8i73UA11A13  -  o.2O3909i65UA11A15  - 
0.0150125181%.^^  +  o.O741071A2uai;la33  -  o.oo8i3ia89UA11A51  +  o.515713852%32  + 
1.726043683%%^  +  0.058004535%%^^  +  o.1744ooo76%3a33  +  o.o6o93o6ooua13a51  + 
1. 217327709%^  -  0. 035947203%%!  +  o.235754i37%5A33  -  0.021220641%%!  - 

0.002222139UA312  +  0 . 0074980141%-% 3  -  0 . 002090749UA31A51  +  O.Ol44l9153UA332  + 

2  3  2 

0.001269098%%!  -  0.000509384%!  -  0.001397249A11  +  o.098283496a11  a13  + 
0.013174031A1;L2A15  +  0.004059203Al12a31  +  0.019404931A-U2A33  +  0.0025l860lA112A51  + 
o.1927118o3A11a132  +  o.725826526a11a13a15  +  o.o496o6oi4a11a13a31  + 
0.076052952A11a13a33  +  o.021958499Al1a13a51  +  0.493373000AuaL52  + 
0.026235560A11A15A31  +  0.1366i4292A11A15A33  +  0 . 017553086a11A15a51  + 

o.oo23858o6a11a312  +  o.oio98318%1a3]a33  +  0.003005606a11a31a51  + 

0  2  3 

o.oo8oio765A11a33  +  o.oo5633548a11a33a51  +  o.ooo98738oa11a51  +  0.603793974a13  - 

1.500490014A132A15  +  o.o43479663A132a31  +  0.324078519A132A33  +  0.02008l073A132A51  + 

4.344o47355A13A152  +  0.i685421i6a13a15a31  +  0. 57978967%%%^  + 

o.o4o344oiia13a15a51  +  o.oo8o85122A13a312  +  0 .019374510A13a31a33  + 

0.009037828a13a31a51  +  o.o6iio94o6a13a332  +  o.oo9834196a13a33a51  + 

0.002764734A13A512  +  i.932393855A153  -  o.260365093A152A31  +  o.8o6221543A152A33  + 

0.060527562A152A51  +  0.007690918a15A312  +  0.03528258oa15a31A33  + 

'2 

o.oioo744i5A15A31A51  +  0.056866015A15a33  +  0.01894594%%%!  + 
0.003407649A15A512  +  0.00044ll37A313  +  0.002235226A312A33  +  0.000853231A31  2v  + 
o.oo2537289A31a332  +  o.oo2683362A31a33a51  +  0.0005698^2%.^  +  o.oo4ii5250A333  + 

0.00i435584A332A51  +  0.000864839A33A512  +  0.000129993A513 
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K5  =  0.009241431U3  -  0.o496o7758ua112  -  1. 70845 3568uAnA13  +  8.125593952UA11A15  - 

0.000252i(62UA11A31  -  0 . 203909l65UA11A33  +  0.0009786001^.^.,^  +  5.9286l99l8UA132  + 
15 . i«525468U2UA13A15  -  0.3^98o8o85UA13A31  +  i.726o43683UA13a33  - 
0.lii-6059989UA13A51  +  10.786l66o48UA152  +  0.855747728UA15A31  +  2.434655418^^33  + 
0.190818970UA15A51  +  0.010124251UA312  -  0.035947203^.^3  +  0.007519578ua31A51  + 
0.ii787705toA332  -  o.02122o64iua33a51  +  0.003141+9810^^  +  0. 0055626974^  - 
0.026496392A112a13  +  i.590352068a112a15  +  o.oi379166oa112a31  +  0.013174031A112a33  + 

0.oo8398672A112a51  +  l.998793077A1]A132  +  5.i29487l40A11Al3A15  + 

0.093686i39A11A13A31  +  o.725826526a11a13a33  +  o.g66i65622A11a13a51  + 

2 

3.547607577A11a15  +  0.7092562i4a11a15a31  +  o.986746oooa11a15a33  + 

0.317769036aua51a51  +  0.006593000A11A312  +  o.02623556oa11a31a33  + 

2 

0 . 00809671  7A11a31a^1  +  o.o68307i46a11a33  +  0 . 017553086ai:_a3  3a?1  + 

0.002597463A11A512  +  2.9o654638oa133  +  25 . 4i42i8o49A132A15  +  o.46oii4o85A132a3;l  - 

2  2*2 
1. 50049001  4a13  a33  +  0.227552472A13  A^  +  32.6l438l245A13A15  + 

i.oo917o478aL3a15a31  +  8.688094670A13a15a33  +  o.42i473114a13a15a51  + 

0 .034424058A  A  ®  +  0.168542116a  A  a  +  0.045551461A  a  a  + 

13  31  13  31  33  13  3i  51 

0.289894836A  ,a  2  +  o.o4o344oiia  a  a_  +  0.015645650A,  a  2  + 

■o  33  13  33  Pi  13  51 

36.413233786aL53  +  o.691341588aL52A31  +  5.797181565^5^33  +  0.287053533^5^  + 
0.120600559A15A312  -  0.520730154a15a31a33  +  o.13468o772a15a31a51  + 
o.8o6221543A15a332  +  o.121055124a15a33a51,  +  o.o4i784273A15A512  +  0.000902957A313  + 
0.007690918a312a33  +  0.00l64888lA312A51  +  o.oi764l290A31A332  + 

0.0100744^2^2^!  +  0.001024491A31A512  +  0. 018955 338a333  +  0.00947297^23^!  + 
0.003407649A33*512  +  o.ooo2i4347A513 
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o  2 

Kg  =  0.003281673U  +  0. 005843 522UA1;l  -  o.oio547763UA11a13  +  o.ooo9786ooua11a15  + 

0.003802360UA11A31  -  0.00813418^^33  -  o.ooi498222UA11a51  +  0.o45884973UA132  - 

o.o46o59989UA13a15  +  0 . 0023680921^  +  0.0609306000^^.^  +  o.0O38i8433UA13a51  + 

0.095409485UA152  +  0.007519578(^^1^  -  0 . 0212206411^^^  +  0.0062899610^^  - 

o.ooo28349.4ua312  -  0.0020907491^3^33  -  0.002157882UA31a51  +  0.0006345490A332  - 

o.ooioi8767UA33a51  -  0.000752348ua512  +  0.003422785A11  +  o.oo89o8319A112a13  + 

o.oo8398672All2a15  +  o.00386o974a112a31  +  0.00251860^^  +  0.002i47644a112a51  4 

0.054554840A  A  2  4  0.066165622A  A  A  +  0.OO8865088A  A  A  4 
11  13’  11  13  15  11  13  3-L 

0.021958499A11AL3A33  4  0. 005494580 A11A13A51  4  o.158884518a1iA152  4 
o.oo8o96717A11aL5a31  4  0 .017553086a11a15A33  4  o.oo5194926auaL5a51  4 
0 . 00l8l66l0A1;LA312  4  0.003005d06A11A31A33  4  0.002213960A11A31A51  4 
0.0028i6774a11a332  4  o.oo197476oal1a33a51  4  0.0007U294a11a512  4  0.028666162^^  4 
0. 227552 472A132A15  4  0.02112808lA132A31  4  0.02008l073A132A33  4  0.012506092A132A51  4 

0.210736557A13AL52  +  0.04555146^^^^  4  0.04034401^^^22  4 

0. 03129130^2 +  0.001927352AL3A312  4  0.009037828A13A31A33  + 
0.002364316a13a31a51  4  o.oo4917098aL3a332  +  0.005^29468aL3a33a51  4 

0.00074170^^^  4  o.095684511A15:>  4  o.o6734o386aL52a31  4  o.o60527562A152A33  4 
o.o4i784273A152A51  4  0.001648881^^2^  +  o.oioo744i5A15a3]a33  + 

0.002048982A15A31A51  4  0.009472971A15A222  4  o.oo6815299A15a33a5i  4 

2  3  2  .  2 

0.00064304^2^5!  +  O.OOO342697A32.  +^0.000853231^3!  A33  +  O.OOO672957A3!  A5!  + 

2  2 
0.001341681^3^33  +  0.00ll39630A3!A33A5!  +  o.ooo45936oA3!A5!  + 

32  23 

0.000478528a33  +  O.OOO864839A33  A^!  +  O.OOO38998OA33A5!  +  0.OOOIO8272A5! 
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TABLE  I 

PARAMETERS  Aij/ao  FOR  FLOW  PAST  A  CIRCULAR  CYLINDER 
AT  VARIOUS  UNDISTURBED-STREAM  MACH  NUMBERS 

&-*] 


A, ,  /a 
11/  0 

A.  _  /a 

13/  0 

A3l/ao 

A._/a 

33/  0 

0.001031 

.001035 

.001103 

.001103 

.001103 

.001105 

-0.0002557 

-.0002557 

-.0002541 

-.0002542 

-.0002541 

-0.0004301 

-.0004301 

-.0004301 

-.0004378 

-0.00001330 

-.00001331 

-.00001337 

0.008605 

.0087^7 

.0093^7 

.009350 

.009346 

.009568 

-0.002199 

-.002198 

-.002140 

-.002141 

-.002135 

-0.003790 

-.003780 

-.003780 

-.004123 

-0.0004862 

-.0004979 

-.0004981 

0.03114 

.03278 

.03520 

.03518 

.03519 

.03540 

-0.008506 

-.008503 

-.007949 

-.007973 

-.007972 

-0.01514 

-.01491 

-.01489 

-.01803 

-0.004733 

-.005045 

-.005051 

0.08307 

.09430 

.1026 

.1018 

.1025 

.1038 

-0.02580 
-.02581 
- .02225 
-.02265 
-.02264 

-0.04693 

-.04381 

-.04297 

-.06157 

-0.03156 

-.03722 

-.03727 

0.1979 

.3242 

.3541 

.3431 

-0.09897 

-.09980 

-.06491 

-0.1657 

.05611 

-0.4292 

0.003155 

.003156  0.03036 
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TABLE  II 

MAXIMUM  VELOCITIES  q/U  OVER  A  CIRCULAR  CYLINDER  COMPUTED  BY 
VARIATIONAL  AND  RAYLEIGH- JANZEN  METHODS  AT  VARIOUS 
UNDISTURBED-STREAM  MACH  NUMBERS 


\  * 

\u 

Variational  method 

Rayleigh-Janzen 
method  (third 
approximation ) 

7 

=  2 

One 

term 

Two 

terms 

. ..  . . - . . .  . 

Three 

terms 

Four 

terms 

Five 

terms 

Si* 

terms 

One 

term 

7  =  1.405 

7  =  2 

0.1 

2.0069 

2.0120 

2.0119 

2.0119 

2.0120 

2.0120 

2.0069 

2.0119 

2.0120 

.2 

2.0287 

2.05H5 

2.0506 

2.0510 

2.0518 

2.0524 

2.0284 

2.0513 

2.0520 

.3 

2.0692 

2.1296 

2.1282 

2.1309 

2.1363. 

2.1364 

2.0681 

2.1314 

2.1360 

.4 

2.1385 

2.2862 

2.2844 

2.2979 

2.3333 

2.3336 

2 . 1040 

2.2836 

2.3027 

.5 

a2.2639 

a2.828l 

a2.827l 
[-  . .  . 

a3.0754 

2.2494 

2.5707 

2.6312 
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TABLE  III 

PARAMETERS  FOR  FLOW  PAS1  A  KAPLAN  BUMP  AT  VARIOUS 

UNDISTURBED -STREAM  MACH  NUMBERS 

E-] 


Mo 

All/ao 

A13/ao 

A3l/ao 

A33/^ 

a15  /ao 

JT 

H 

O 

0.5 

0.009459 

-0.003325 

-0.001799 

0.000558 

0.000180 

0.000156 

0.75 

0.05289 

-0.02119 

-0.01607 

0.01662 

0.02490 

-0.01327 

0.83 

0.10501 

-0.04533 

-0.02624 

0.04201 

0.00799 

-0.07141 

0.90 

0.092073 

.20799 

.2554 

.2736 

-0.081943 

-.08626 

-.1064 

-0.1917 

-.2256 

0.1558 

VELOCITY  DISTRIBUTION  qjU  OVER  SURFACE  OF  A  KAPLAN  BUMP  COMPUTED  BY  VARIATIONAL  AND 
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Four  parameters  of  the  series  taken. 
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TABLE  VI 

PRESSURE  COEFFICIENTS  COMPUTED  BY  VARIATIONAL  METHOD 
FOR  VARIOUS  UNDISTURBED-STREAM  MACH  NUMBERS  AT 
VARIOUS  VALUES.  OF  X  ON  BUMP 


Vp 

x  n'nv 

O 

It 

O 

as 

Mq  =  0.5 

. 

. 

M0  =  0.75 

Mq  =  0.83 

0 

-0.168777 

-0.196976 

-0.278837 

-0.380572 

.1 

-.164025 

-.191428 

-.270078 

-.366493 

.2 

-.151114 

-.176112 

-.245300 

-.327481 

.3 

-.130820 

-.151796 

-.206697 

-.266558 

.4 

-.101030 

-.117012 

-.148537 

-.187512 

.5 

-.067709 

-.077888 

-.095925 

-.102108 

.6 

-.030022 

-.033808 

-.033924 

-.017223 

•  7 

.010174 

.012772 

.027712 

.058712 

.8 

.051714 

.060936 

.085342 

.118537 

.9 

.093315 

.107692 

.135733 

.155920 

.975 

.124466 

.143040 

.167266 

.161538 
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TABLE  VII 

PRESSURE  COEFFICIENTS  COMPUTED  BY  VARIOUS  METHODS 
FOR  Mq  =  0.83  AT  VARIOUS  VALUES  OF  X  ON  BUMP 


\Cp 

X  >v 

Pr  andt 1 -Glauert 
method 

Von  Karman-Tsien 
method 

Perturbation 

method 

Variational 

method 

0 

-0.30245 

-0.324138 

-0.35564 

-0.380572 

.1 

-.291*033 

-.314480 

-.34310 

-.366493 

.2 

-.270905 

-.288167 

-.30771 

-.327481 

.3 

-.234545 

-.247375 

-.25271 

-.266558 

.k 

-.181081 

-.188634 

-.18424 

-.187512 

.5 

-.121378 

-.124726 

-.10861 

-.102108 

.6 

-.053787 

-.054434 

-.03270 

-.017223 

.7 

.018234 

.018160 

.03827 

.058712 

.8 

.092715 

.090853 

.10222 

.118537 

•  9 

.167303 

.161335 

.16275 

.155920 

.975 

.223151 

.212658 

.21232 

.161538 

Velocity  at  surface  of  cylinder,  q/U 
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0  10  20  30  40  50  60  70  80  90 


Angle,  0,  deg 

Figure  2.-  Velocity  distribution  over  a  circular  cylinder  at  Mq  =  0.4 

computed  by  variational,  Rayleigh-Janzen,  and  incompressible -flow 
methods. 
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Maximum  local  Mach  number 
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0  .1  .2  .3  A  .5 

Undisturbed-stream  Mach  number,  MQ 

Figure  3.-  Maximum  local  Mach  number  at  various  undisturbed-stream  Mach 
numbers  computed  by  variational  and  Rayleigh-Janzen  methods. 


Local  Mach  number  at  surface  of  cylinder 
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Angle,  e ,  deg 

Figure  4.-  Local  Mach  number  over  surface  of  cylinder  computed  by 
variational  and  Rayleigh- Janzen  methods  for  various  undisturbed- 
stream  Mach  numbers. 
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Method 


